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FOREWORD 

This  report  presents  the  results  of  work  which  was  supported  by  the  Naval  Air 
Systems  Command,  “Advanced  Reliable  Solid  Propulsion  Design,”  Propulsion 
Administrators,  Mr.  Robert  Heitkotter  (NAVAIR)  and  Dr.  Charles  Thelen  (NWC)  and 
by  the  Naval  Sea  Systems  Command,  “Plume  Low  Signature  Requirements  for 
Point/Area  Defense  Solid  Propellant  Motors,”  J.  Murrin  and  R.  Cassel.  The  Program 
Manager  was  Mr.  George  Buckle,  Naval  Ordnance  Station,  Indian  Head,  MD. 

This  report  is  a  revision  of  NWC  TM  3638.  Since  the  original  report  was 
published  in  November  1978,  an  error  in  the  computer  program  that  was  used  to 
generate  much  of  the  data  in  the  report  has  been  corrected.  As  a  result.  Tables  2 
throught  15  and  the  computer  program  in  the  Appendix  have  been  changed.  Also,  a 
slight  change  has  been  made  in  the  percent  value  in  Figure  4. 

This  report  is  issued  for  use  at  the  working  level  only  and  does  not  reflect  an 
official  position  of  the  Naval  Weapons  Center. 


Duane  H.  Williams 
Head,  Advanced  Technology  Division 
Ordnance  Systems  Department 
28  January  1979 


NWC  TM  3638,  Revision  1,  published  by  Code  324,  30  copies. 
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ABSTRACT 

Under  certain  climatic  conditions  rocket  motor  exhaust  plumes  can  form 
contrails  which  are  highly  visible  at  long  ranges.  It  is  difficult  to  predict  accurately  the 
formation  and  visibility  of  these  contrails  unless  the  important  climatic  factors 
(temperature  and  moisture  content)  are  known.  This  report  proposes  a  preliminary 
definition  of  a  climate  at  40°N  and  50°  N  latitude  for  altitudes  from  2-15  km  based 
on  a  review  of  relevant  literature.  Both  temperature  and  moisture  content  are 
expressed  in  terms  of  mean  values  and  standard  deviations  so  that  the  probability  of 
occurrence  of  any  climatic  conditions  can  be  estimated. 
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INTRODUCTION 

Moisture  laden  exhaust  contrails  can  greatly  increase  the  detectability  of  missiles 
and  aircraft.  The  problem  is  most  severe  for  those  rocket  motors  which  contain,  as 
products  of  combustion,  water  soluble  gas  and  particle  species  which  increase  the 
propensity  for  moisture  condensation.  Therefore,  avoidance  of  contrail  formation  is  an 
important  design  criterion  for  certain  missions  required  of  modem  missile  propulsion 
systems.  Changes  made  in  missile  propellant  ingredients  in  order  to  meet  this  criterion 
often  result  in  reduced  rocket  motor  total  energy.  To  minimize  this  compromise  of 
motor  energy  it  is  important,  in  the  system  design  phase,  to  perform  trade-off  analyses 
of  the  effects  of  energy  and  contrail  formation  on  the  performance  of  a  missile’s 
mission. 

Previous  studies  have  yielded  techniques  for  predicting  missile  performance1  and 
the  visibility  of  primary  smoke*  trails.  Although  techniques  have  also  been  developed 
foi  predicting  the  formation  and  visibility  of  secondary  smoke**  contrails3,4,  these 
predictions  require  inclusion  of  ambient  temperature  and  humidity  values. 

Temperature  and  humidity  values  needed  for  incorporation  in  contrail 
predictions  may  be  obtained  in  several  ways: 

1.  Use  specific  values  based  on  specific  measurements 

2.  Use  standard  climatic  values5,6 


*  “Primary  smoke"  is  composed  of  the  solid  effluent  from  propulsion  systems. 

**  “Secondary  smoke"  is  formed  by  condensation  of  plume  and  atmospheric  water  on  the  microscopic 
particle  efflux;  it  often  contains  water  soluble  species  too. 

*  Naval  Weapons  Center.  Tactical  Missile  Encounter  (TAME)  Simulation  Computer  Program,  by  D.  R. 
Peters,  R.  D.  Feldman,  H.  Sobel  and  M.  J.  Stevens.  China  Lake,  Calif.,  NWC,  March  1970.  (NWC  TP  4866, 
publication  UNCLASSIFIED.) 

*  A.  C.  Victor  and  S.  H.  BreU.  “A  Simple  Method  for  Predicting  Rocket  Exhaust  Smoke  Visibility,”  J 
SPACECRAFT  AND  ROCKETS ,  Vol.  14,  No.  9  (September  1977),  pp.  52M33. 

3  Naval  Weapons  Center.  "Computer  Codes  for  Predicting  Secondary  Smoke  Formation  in  Free  Jets  and 
Smoke  Chambers,"  by  A.  C.  Victor.  China  Lake,  Calif.,  NWC,  February  1978.  (NWC  TM  3361,  publication 
UNCLASSIFIED.) 

*  H.  Hoshizaki,  et.  al.  “Plume  Visibility  Detection  Study,”  presented  at  the  10th  JANNAF  Plume 
Technology  Meeting,  1417  September  1977.  (NWC  TM  3361,  publication  UNCLASSIFIED.) 

3  U.S.  Standard  Atmosphere  1962.  U.S.  Government  Printing  Ofllce,  December  1972,  Washington,  D.C. 

6  Air  Force  Cambridge  Research  Laboratories,.  Handbook  of  Geophysics  and  Space  Environments,  S.  L. 
Valley,  Ed.  Bedford  Maw.,  1965,  Chapter  3.  (UNCLASSIFIED.) 
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3.  Use  tabulated  climatic  values  (footnotes  5  and  6)  and  make  the  contrail 
predictions  for  a  world-wide,  year  round,  range  of  climates 

4.  Use  a  statistically  defined  climate  model  (humidity,  temperature,  altitude) 
which  corresponds  to  observations  and  makes  the  contrail  prediction  at 
statistically  significant  points  in  the  climate 

Each  of  these  methods  of  selecting  climatic  values  has  distinct  advantages  and 
disadvantages.  The  first  method  is  useful  only  to  verify  the  predictive  technique  by 
comparison  with  previously  made  measurements;  it  has  no  relevance  for  system  design 
studies. 


The  second  method  is  useful  for  a  quick  look  approach,  but  it  cannot  provide 
probabilities  of  contrail  formation  which  are  needed  in  the  design  study. 

The  third  method  is  capable  of  generating  the  predictions  needed  in  design 
studies,  but  at  an  excessive  cost,  because  of  the  many  climate  conditions  involved.  This 
method  is  most  useful  for  evaluating  contrail  formation  in  specific  operation  theaters. 

The  fourth  method  is  essentially  a  condensation  of  the  third  method  in  that 
weather  statistics  are  included  in  the  climate  model.  Thus  with  only  a  few  calculations 
the  probability  of  contrail  formation  for  a  specific  missile  propulsion  system  can  be 
predicted  as  a  function  of  altitude  and  latitude. 

This  report  presents  a  simple  mode!  climate-  based  on  the  fourth  method  above, 
which  averages  annual  climate  variations  at  specific  latitudes.  Latitudes  of  40°N  and 
50°N  were  chosen  as  being  typical  of  the  United  States  and  Western  Europe, 
respectively. 

A  more  complex  statistical  model,  based  on  accumulated  world-wide  statistics, 
although  desirable,  was  beyond  the  manpower  and  funding  limits  of  the  projects  which 
supported  this  study.  Such  a  model  would  be  capable  of  predicting  weather 
probabilities  on  any  specific  date,  whereas  the  model  described  in  this  report  can  only 
predict  probabilities  for  an  arbitrary  date. 


CLIMATE  MODEL 

The  climate  model  needed  for  contrail  prediction  has  four  elements: 

1.  Mean  temperature  as  a  function  of  altitude 

2.  Standard  deviation  of  temperature  as  a  function  of  altitude 

3.  Mean  humidity  as  a  function  of  altitude 

4.  Standard  deviation  of  humidity  as  a  function  of  altitude 
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TEMPERATURE  MODEL 

Temperature  as  a  function  of  altitude  has  been  measured  carefully  for  many 
years,  and  sufficient  data  are  reported  to  define  both  elements  1  and  2  above.  From 
Table  3-9  in  Handbook  of  Geophysics  and  Space  Environments  (footnote  6),  Equations 
1  and  2  have  been  generated  to  represent  mean  annual  temperature,  in  degrees  Celsius, 
at  40°N  and  50°N  latitude  respectively,  as  functions  of  altitude  (y)  in  km  up  to  20 
km. 


T40oN  =  28.56-  10.618y  +  0.3168y2  (1) 

T50°n  =  12.95  -  8.774y  +  0.2738y2  (2) 

Equations  3  and  4  represent  the  January  and  July  mean  variations  of  temperature 
about  the  annual  mean: 

t40->N  =  8.0-  0.56y  +  0.01 53y2  (3) 

t50°N  =  10  3i  -  °-7y  +  0.0167y2  (4) 


The  mean  temperature  in  January  and  July  are  given  by: 

TGan)40»N  =T4qon  -t40°N 

=  T40°N  +  t40°N 
=  T40°N 

T(july)40°N  =  T40.n  +  t40«N 

=  T40°N  "  t40°N 
=  T40°N 

TG'an)50®N  =T5o°n  “t50°N 

TGuly)50°N  =T50°N  +  t50°N 

The  standard  deviations  (a)  about  the  January  and  July  mean  values  (up  to  15  km) 
are  approximately: 

oGan)40»N  =  ±5°C 
oGuly)40»N  =  ±3°C 

t  °Oan)50°N  =  *6  C 

oGuly)5o°N  =  ±4°C 


0  to  12  km  altitude 
14  to  20  km 
12  to  14  km 

0  to  1 2  km 
14  to  20  km 
12  to  14  km 

0  to  20  km 
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HUMIDITY  MODEL 

The  information  in  footnote  6  and  more  recent  measurements7  show  that  mean 
atmospheric  moisture  content  can  be  well  represented  by  Equation  9  for  the  mixing 
ratio  W  =  mass  of  atmospheric  water/mass  air,  in  parts  per  million, 

W  =  10<3-7  ~  °-192y)  for  0-10  km  altitude 

=  104"1 3  “  °.235y)  for  10-15  km  altitude  (9) 

up  to  an  altitude  of  15  km  (about  50,000  feet).  The  data  in  footnote  7  indicate  that 
the  fractional  standard  deviation  about  the  means  is  about  0.6  at  all  altitudes. 

The  actual  water  vapor  pressure  (ey)  is  obtained  from  the  mixing  ratio  by 
Equation  10 


PWX  10“6  PW 

e  - -  s -  nn'i 

v  0.622 +  WX  10-6  0.622  X  106 

in  the  units  of  atmospheric  pressure  (P). 

Standard  atmospheric  pressure  (footnote  6)  as  a  function  of  altitude  is  given  in 
Table  1. 


TABLE  1.  Atmospheric  Pressure  (P)fl  as  a 
Function  of  Altitude  (y). 


y,  km 

P,  mb 

y,  km 

P,  mb 

0 

1013 

12 

193 

2 

795 

14 

141 

4 

616 

16 

103 

6 

472 

18 

75.0 

B 

356 

20 

54.7 

10 

264 

22 

40.0 

e  The  following  constants  may  be  used  to 
convert  to  other  pressure  units: 

mb  x  0.00098692  =  atmospheres 

mb  x  10®  -  pascals  (N  m-2) 

mb  x  0.014504  -  pounds/inch2 

mb  X  0.7502  =  mm  of  mercury  at  0°C  1 


7  Air  Force  Cambridge  Research  Laboratories.  “Humidity  Up  to  the  Mesopause by  N.  Sissenwine,  D.  D. 
Granthan  and  H.  A.  Sametx  Bedford,  Mass.,  October  1968.  (AFCRL-68-0S  JO,  AO  679996,  UNCLASSIFIED.) 
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By  combining  Equations  9  and  10,  Table  1,  and  water  vapor  pressure  curves 
(footnote  3,  Handbook  of  Chemistry  and  Physics*  or  Figure  2),  one  can  describe  the 
humidity  in  terms  of  a  frost  (or  dew)  point.  This  is  shown  in  Figure  3. 


TEMPERATURE,  *C 


FIGURE  2.  Temperature  Dependence  of  the  Equilibrium  or  Saturated 
Vapor  Pressure  of  Pure  Water. 


8  The  Chemical  Rubber  Co.  Handbook  of  Chemistry  and  Physics.  R.  C.  Wets!  and  S.  M.  Selby,  Ed.,  47th 
ed.  (1966),  pp.  D-105,  106.  ' 
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It  is  interesting  to  note  that  if  one  overlays  a  transparency  of  Figure  3  on 
Figures  la  and  lb,  some  information  on  the  probability  of  cloud  formation  emerges.  If 
the  frost  point  is  greater  than  the  ambient  temperature,  clouds  will  form. 

The  relative  humidity  at  a  temperature  is  the  ratio  of  the  water  vapor  pressure 
to  the  saturation  pressure  of  water  at  that  temperature  (ps(T)).  If  the  calculated 
relative  humidity  is  greater  than  unity,  clouds  will  form.  In  the  presence  of  clouds, 
contrail  formation  does  little  to  increase  the  visibility  (detection  probability)  of  a 
missile.  Therefore  we  are  only  interested  in  those  portions  of  Figures  1  and  3  for 
which  the  frost  point  is  lower  than  the  ambient  temperature. 


FIGURE  3.  Frost  (Dew)  Point  Temperature  vs.  Altitude  from  Humidity 
Model.  (Fractional  standard  deviation  is  ±0.6  of  W  or  ev.) 
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COMPOUND  MODEL 

In  order  to  move  ahead  from  this  stage  of  the  model,  two  assumptions  are 
necessary:  (1)  the  relationship  between  the  temperature  and  humidity  models  and  (2) 
the  shapes  of  the  temperature  and  humidity  probability  distributions. 

Independence  of  the  temperature  and  humidity  is  assumed,  and  the  shapes  of 
the  probability  curves  are  assumed  to  be  “triangular”  as  shown  in  Figure  4.  These 
assumptions  were  made  in  lieu  of  readily  accessible  contradictory  data. 

Following  this  assumption,  it  is  possible  to  calculate  cumulative  probabilities  of 
any  value  of  temperature  and  humidity  and  of  the  joint  probability  of  any 
temperature-humidity  pair.  The  probability  of  real  interest  is  the  joint  probability  that 
temperature  will  be  below  a  specified  value  and  that  humidity  will  be  above  a  certain 
value.  The  joint  probabilities  and  relative  humidities  for  this  set  are  shown  in  Tables  2 
through  15  for  the  altitude  range  of  2-14  km  at  both  40°N  and  50°N  lattitude.  These 
tables  were  generated  by  the  computer  program  given  in  the  Appendix. 

The  number  pairs  enclosed  by  the  rectangles  in  Tables  2  through  15  represent 
an  11  by  11  matrix  of  joint  probabilities  and  relative  humidities.  The  upper  number  in 
each  pair  is  the  joint  probability;  the  lower,  the  relative  humidity.  The  joint 
probability  represents  the  probability  that  temperature  is  lower  than  the  value  in  the 
far  left  column  (outside  the  rectangle)  and  that  humidity  is  greater  than  the  value  in 
the  top  row  (outside  the  rectangle).  The  right-most  column  in  each  table  is  the 
probability  that  the  temperature  is  less  than  the  value  in  the  left-most  column.  The 
lowest  row  in  each  table  is  the  probability  that  the  humidity  (ev  in  mm  Hg)  is  greater 
than  the  value  in  the  top-most  row.  The  second  column  in  each  table  is  the  saturated 
vapor  pressure,  ps(T)  or  PSAT,  associated  with  the  temperature,  TC,  in  the  first 


FIGURE  4.  Assumed  Probability  Curve 
Shape  for  Temperature  and  Humidity. 
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column.  The  relative  humidity  is  just  the  ratio  ey/ps(T).  Regions  in  Tables  2  through 
15  where  the  relative  humidity  is  greater  than  unity  correspond  to  natural  cloud 
formation. 

For  example,  if  one  looks  at  the  central  pair  (the  6,6  matrix  position)  in  the 
matrix  of  Table  2,  the  joint  probability  is  0.25;  the  relative  humidity  is  0.256  (or 
25.6%).  This  pair  corresponds  to  a  temperature  of  7.6°C  (which  has  a  saturated  vapor 
pressure  of  7.76  mm  Hg)  and  an  actual  vapor  pressure  of  ey  =  1.98  mm  Hg.  There  is 
a  0.5  probability  that  temperature  is  less  than  7.6°C  at  this  altitude  and  latitude  and 
also  a  0.5  probability  that  the  humidity  is  greater  than  1.98  mm  Hg. 


CONCLUSIONS 

A  simple  model  for  defining  the  atmospheric  temperature  and  humidity  at 
regions  of  operational  interest  for  tactical  missiles  in  the  temperate  zone  of  the  earth’s 
northern  hemisphere  has  been  described  in  this  report.  With  this  model  one  can  define 
a  probability  distribution  for  temperature  and  humidity  on  any  arbitrary  day. 

Definition  of  probabilities  of  these  atmospheric  variables  may  also  be  desirable 
for  any  particular  day.  In  order  to  obtain  appropriate  values  for  this  approach  it  will 
be  necessary  to  computerize  extensive  records  of  atmospheric  fluctuations  and  derive 
appropriate  statistical  and  perhaps  analytical  approaches  which  are  beyond  the  funding 
limits  of  the  projects  which  supported  this  study. 


n 


TABLE  2.  Climate  for  40°N  Latitude 
2.0  km  Altitude  795.0  MB  Pressure 
(Joint  Probability /Relative  Humidity). 

EV,  MR  HG 
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862  .807  .7 *2  .670  .569  .500  .(11  .350  .258  .194  .138 

PROBABILITY  THAT  HUH10ITY  .GT.  EV 


TABLE  3.  Climate  for  40°N  Latitude 
4.0  km  Altitude  616.0  MB  Pressure 
(Joint  Probability/Relative  Humidity). 
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862  .807  .742  .670  .569  .530  .411  .330  .255  .194  .138 

PROBABILITY  that  NURIDITY  .ST.  EV 
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TABLE  13.  Climate  for  50° N  Latitude 
10.0  km  Altitude  264.0  MB  Pressure 
(Joint  Probability /Relative  Humidity). 
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TABLE  15.  Climate  for  50° N  Latitude 
14.0  km  Altitude  141.0  MB  Pressure 
(Joint  Probability /Relative  Humidity). 
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Appendix 

COMPUTER  PROGRAM  FOR  GENERATING  PROBABILITY  TABLES 
c 

C**»  PR  06  RAM  TO  CALCULATE  CLIMATE  AT  40  AND  50  DEG  NORTH  LATITUDE 
C***«CORRECTED  VERSION,  4  JAN  1979 
C 

C  Y  IS  ALTITUDE,  KM 

C  X  LAT  IS  LATITUDE,  DE6  N 

C  PATH  IS  ATMOSPHERIC  PRESSURE  IN  MB 

C 

DIMENSION  H(11),PH(11),PT(11),PJ(11),RH(11),PHH(11),PTT(11) 
5  READC5, 100, END-75)  T,XLAT,PATM 
WRITE (6,200) XLAT,Y,P ATM 
IFtXLAT.LT. 50.)  toO  TO  .10 
C 

C  FOR  LAT  *  SO  DEG  N 

C 

TEMP=12.95-8.774»Y*.2738*Y**2. 

TE  *  1Q.3-0.7*Y*.0167*Y**2. 

SI6U  *  6. 

S16L  *  4. 

TU*  TEMP*TE*SI6L 
TL  *  TEMP-TE-SI6U 
GO  TO  20 

C  FOR  LATITUDE  =  40  DEG  N 

10  TEMP  =  28. 56-10. 618*Y*.3168*Y**2. 

TE  *  8.  -  . 5d*y  ♦ .0153*Y**2 . 

SIGU  -  5. 

S1GL  =  3. 

IF(Y.GT.14.)  60  TO  IS 
TU  *  TEMP  ♦  TE  ♦  SI GL 
TL  *  TEMP  -  TE  -  SIGU 
GO  TO  2C 

15  TU  *  TEMP  ♦  TE  ♦  SIGU 
TL  *  TEMP  -  TE  -  SIGL 
C  CALCULATE  HUMIDITY 

20  IFtY.faT.IO.)  GO  TO  21 

M  *  10. **(3.7  -  . 192 *Y) 

GO  TO  22 

21  N  »  10. **(4.13  -  ,235*Y) 

22  EV  *  PATH  *  M/622000. 

EVH6  =  EV*  .7502 
EVHGU  *  1 • 6*E VHG 
EVH6L  *  .4*EVHG 

C  CALCULATE  TEMPERATURE  STATISTICS 

BT  *  (TU  -  TD/.72 
AT  =  TL  -  .  1 4 *BT 

C  CALCULATE  HUMIDITY  STATISTICS 

BH  *  -(EVHGU  -  EVMGD/.72 
AH  *  EVHGU  -  .  14*BH 

C  CALCULATE  JOINT  PROBABILITIES  AND  RELATIVE  HUMIDITIES 

PH( 1 )  *  .8625 
PHM ( 1 )  *  .86 
DO  30  1  *  1,11 
H (I )  *  AH  ♦  PHH(I)*BH 
I  Fd.LE.5)  DED  * (0.8 25* ( 1-1 )  *5.6>*.01 
IFU.&T.S)  DED  «(0.825*(10-I)  *  5.6)*. 01 
PH ( I* 1)  *  PM( I )  -  DED 
PHH(1*1)  *  PHH(I)  -  .072 
30  CONTINUE 

WRITE (6,2 1 0) 

WRITE (6,220)  (H(I),I*1,11) 
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61*  PT(1)  *  .8625 

62*  PTT(1>  =  .86 

63*  00  40  J*1,11 

64*  TC  '  AT  ♦  PTT  C  J )  * BT 

65*  T  =  TC  ♦  273. 

66*  IFCT.LT. 256.)  60  TO  520 

67*  510  VPS  =  (  ( (  .  1 35601  i-?*T-.1 5 1036E-4)  »T*.00618156>»T-1 .07005)*  T+64  .058 

68*  60  TO  530 

69*  520  VPS*( (( .538009E-05-.473762E  -03*T >*T- .00 243 188 ) *T* . 54582?) *T-S0. 27 

7C*  189 

71*  530  VPS  =  10.»*VPS 


72* 

73* 

74* 

75* 

76* 

77* 

78* 

79* 

80* 

81* 

82* 

83* 

84* 

85* 

86* 

87* 

88* 

89* 

90* 

91* 

92* 

53* 

94* 

95* 

96* 

97* 

98* 


PH(1>  *  .8625 
00  35  I  *1,11 
PJCI)  *  PHCI)*PTCJ> 

RH(I)  *  HCD/VPS 
35  CONTINUE 

WRITE(6,230)TC,(PJ(I),X*1 ,11) ,PT(J> 

WRITE(6,240)VPS,CRH(I),I*1,11> 

1FCJ.LE.5)  0£0  =(0.825«(J-1>  *5. 6)*. 01 
IFU.6T.5)  0E0  =C0.825*C10-J)  ♦  5. 6)*. 01 
PTCJ+1)  *  PTCJ)  -  DEO 
PTTCJ  +  1)  *  PTT(J)  -  .072 
40  CONTINUE 

WRITE (6 ,250)  (PH ( I) , 1*1 ,11) 

WRITE (6,260) 

6  0  TO  5 
75  CALL  EXIT 
100  FORMATO 

200  f  ORH  AT(  1H1  T4 1  , 'CLIMATE  FOR  ',F4.1,'DEG.  N  LATITUDE'  /T40,F  4.1,  'KN 
1ALTITU0E', 56.1 ,'  MB  PRESSURE',  /T40, '(JOINT  PROBABILITY  /  RELATIVE 
2  HUMIDITY  )'/T45, '************* *•****•'/ T50,'€V,  HN  HG') 

210  for«at<iho,'tc',t9,'psat',ti09,'prob  temp  is') 

220  FORMATClH  ,'OtG  C',T9,'NH  H6' ,T17 .1 1 E8 .3  ,6X 'LT  TC') 

230  FORMAT  (  / 1 M  F6 . 1 ,T 1 7, 1 1 E8. 3 , F8 .3 > 

240  FORMATClH  T9.1ZEE.3) 

250  F0RMAT(1H0T15,11F8.3) 

260  FORHATdH  T40, 'PROBABILITY  THAT  HUMIDITY  .61.  EV'> 

END 
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